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Autocorrelation-Based  Spectrum  Sensing 
Algorithms  for  Cognitive  Radios 

Takeshi  Ikuma  and  Mort  Naraghi-Pour 
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Louisiana  State  University 
Baton  Rouge,  LA  70803 
Email:  {tikuma@lsu.edu,  mort@ece.lsu.edu} 


Abstract — Cognitive  radio  is  an  enabling  technology  for  op¬ 
portunistic  spectrum  access.  Spectrum  sensing  is  a  key  feature 
of  a  cognitive  radio  whereby  a  secondary  user  can  identify 
and  utilize  the  spectrum  that  remains  unused  by  the  licensed 
(primary)  users.  Among  the  recently  proposed  algorithms  the 
covariance-based  method  of  [1]  is  a  constant  false  alarm  rate 
(CFAR)  detector  with  a  fairly  low  computational  complexity. 
The  low  computational  complexity  reduces  the  detection  time  and 
improves  the  radio  agility.  In  this  paper,  we  present  a  framework 
to  analyze  the  performance  of  this  covariance-based  method.  We 
also  propose  a  new  spectrum  sensing  technique  based  on  the 
sample  autocorrelation  of  the  received  signal.  The  performance  of 
this  algorithm  is  also  evaluated  through  analysis  and  simulation. 
The  results  obtained  from  simulation  and  analysis  are  very  close 
and  verify  the  accuracy  of  the  approximation  assumptions  in 
our  analysis.  Furthermore,  our  results  show  that  our  proposed 
algorithm  outperforms  the  algorithm  in  [1]. 

Index  Terms — Spectrum  sensing.  Dynamic  Spectrum  Access, 
Cognitive  Radio,  Autocorrelation-based 

L  Introduction 

In  2002,  Federal  Communications  Commission  (FCC)  Spec¬ 
tral  Policy  Task  Force  reported  that  typical  radio  channel 
occupancy  is  less  than  15%  while  the  peak  occupancy  is 
close  to  85%  [2].  To  increase  the  spectrum  usage  efficiency, 
the  task  force  recommended  the  development  of  opportunistic 
or  dynamically  spectrum  access.  Dynamic  spectrum  access 
meshes  well  with  the  cognitive  radio  [3]  paradigm  which, 
as  an  evolution  of  software-defined  radios,  is  aware  of  its 
surrounding  environment  and  can  accordingly  adapt  its  internal 
operating  states  [4]. 

Dynamic  spectrum  access  requires  frequency  agile  radios 
that  can  monitor  and  identify  the  spectral  bands  that  are  unused 
by  the  primary  licensed  users  (the  so-called  white  spaces). 
The  radio  will  then  dynamically  adjust  its  carrier  frequency, 
transmit  power,  modulation,  coding,  etc,  in  order  to  make  best 
use  of  the  available  spectrum  and  achieve  the  desired  quality 
of  service.  In  order  to  avoid  interference  to  primary  users, 
spectrum  sensing  must  detect  signals  with  low  probability  of 
detection.  Moreover,  low  computational  complexity  and  ease 
of  implementation  is  required  in  order  to  facilitate  radio  agility. 

Akyildiz  et  al.  categorized  non-cooperative  spectrum  sens¬ 
ing  into  three  categories  [5]:  energy  detection,  matched-filter 
detection,  and  cyclostationary  feature  detection.  Spectrum 
sensing  using  energy  detection  is  easy  to  implement  and 


performs  well  if  the  noise  power  at  the  receiver  is  known. 
However,  uncertainty  in  the  noise  power  can  significantly 
degrade  the  performance  of  energy  detectors.  Furthermore, 
estimation  of  the  noise  power  which  is  required  in  such 
cases  leads  to  the  so  called  “SNR  wall”  phenomena  [6]. 
The  detectors  based  on  matched  filtering  and  cyclostationary 
features,  on  the  other  hand,  rely  on  the  a  priori  knowledge 
of  the  signal  parameters.  The  matched-filter  based  system  is 
a  coherent  system  and  must  have  a  priori  knowledge  of  the 
modulation  type  and  carrier  frequency  of  the  primary  user.  The 
cyclostationary  feature  detectors  also  require  some  parameters 
of  the  primary  signal  such  as  symbol  rate  and  have  high 
computational  requirements. 

An  alternative  approach  is  autocorrelation-based  method 
suggested  in  [1].  Their  approach  relies  on  the  fact  that  when 
the  receiver  bandwidth  is  greater  than  the  bandwidth  of  the 
signal,  the  autocorrelation  function  of  signal  plus  noise  is 
distinguishable  from  that  of  noise  alone  (which  is  assumed 
to  be  white).  Because  autocorrelation  function  is  a  one¬ 
dimensional  entity  as  opposed  to  the  two-dimensional  cyclic 
autocorrelation  function,  the  burden  in  computing  the  decision 
statistics  is  drastically  reduced.  The  decision  statistic  in  [1] 
is  based  on  the  Frobenious  norm  of  the  correlation  matrix 
of  the  received  signal.  They,  however,  have  not  analyzed  the 
performance  of  their  detector  and  only  presented  an  empirical 
false  alarm  rate. 

In  this  paper,  we  present  an  accurate  analysis  of  the 
autocorrelation-based  detector  in  [1]  using  the  results  in  [7] 
on  the  first-order  autoregressive  (AR)  spectral  estimator.  The 
presented  analysis  can  be  extended  for  the  performance  of 
other  autocorrelation-coefficient  based  algorithms.  We  also 
propose  a  new  autocorrelation-based  detector  by  considering 
additional  properties  of  the  autocorrelation  function  of  typical 
baseband  communication  signals.  These  properties  greatly 
simplify  the  complexity  of  the  detection  statistics  resulting  in 
a  linear  detector  as  opposed  to  the  quadratic  form  in  [1].  We 
also  evaluate  the  performance  of  our  detector  through  analysis 
and  simulation. 

The  remainder  of  this  paper  is  organized  as  follows.  In 
Section  II  we  present  the  spectrum  sensing  problem,  the  sys¬ 
tem  configuration  and  the  assumptions  on  the  received  signal. 
Section  III  introduces  the  detector  in  [1]  as  well  as  our  new 
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r{t) 
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sampling  system  . 

Fig.  1.  Simplified  block  diagram  of  ideal  RF  front-end. 


detector.  Section  IV  presents  the  analysis  of  the  performance 
of  detectors  in  terms  of  detection  and  false  alarm  probabilities. 
Section  V  presents  numerical  results  from  simulation  and 
analysis.  Finally,  concluding  remarks  are  provided  in  Section 
VII. 

II.  Spectrum  Sensing  System 

For  the  purpose  of  spectrum  sensing  the  cognitive  radio 
front-end  can  be  simplified  as  shown  as  in  Fig.  1.  The  radio 
receives  an  RF  signal  r(t)  and  after  down  conversion,  low-pass 
filtering,  and  sampling,  obtains  the  complex  baseband  signal 
{xn}-  The  target  center  frequency  is  Hz,  the  bandwidth  of 
the  low-pass  filter  is  {— fbw,  fbw)  Hz,  and  the  sampling  rate 
Ts  is  given  by  =  {2fbw)~^- 

The  complex  baseband  signal  x{t)  is  modeled  as  x{t)  = 
+  n{t)  where  s{t)  is  the  primary  baseband 
communication  signal,  n{t)  is  the  complex  noise  process, 
and  where  and  9q  denote  the  frequency  and  phase  offsets 
of  the  local  oscillator  from  that  of  the  primary  transmitter, 
respectively.  The  value  of  77  €  {0, 1}  determines  the  presence 
or  absence  of  the  primary  signal  s{t).  Therefore  the  detection 
of  the  primary  signal  is  described  by  the  following  binary 
hypotheses  testing  problem. 


Hq:  77  =  0, 

i^i:  17  =  1, 


primary  signal  absent 
primary  signal  present 


(1) 


(2) 


signal-to-noise  ratio  (SNR)  of  a;„  is  denoted  by 

A  rs,o 

3'  =  ~- 

Assuming  that  {s„}  and  {77„}  are  uncorrelated,  we  express 
the  conditional  autocorrelation  function  of  as 


(4) 


It  is  further  assume  that  the  real  and  imaginary  parts  of  {s„}, 
namely  {5ft(s„)}  and  {Sr(s„)}  are  independent.  The  following 
lemma  results. 

Lemma  1:  If  the  real  and  imaginary  components  of  a 
complex-valued  zero-mean  WSS  random  process  are  indepen¬ 
dent,  then  the  autocorrelation  function  of  the  process  is  an 
even  function. 

Proof:  Let  a  complex-valued  zero-mean  WSS  process 
Xn  =  Un  +  jzn  where  T/n  and  Zn  are  mutually  independent 
real  WSS  processes.  Then, 

Xxx^l  —  P[XnXn—l\ 


=  E[{yn  +  jZn){yn-l  “  j Zn-l)] 


(5) 


where  *  denotes  complex  conjugation.  Since  a;„  is  zero-mean 
and  7/„  and  are  independent,  we  have 


rxx,l  =  E[ynyn-l]E[ZnZl_i] 


(6) 


The  primary  signal  s{f)  is  unknown  and  is  modeled  as  a 
complex-valued  zero-mean  wide-sense  stationary  (WSS)  pro¬ 
cess,  characterized  by  its  autocorrelation  function  rs{T)  = 
i?[s(f)s(f  —  r)*].  Furthermore,  s{t)  is  band-limited  in  the 
frequency  range  {-fb,fb)  Hz  where  fb  <  fbw 

The  spectrum  sensing  algorithm  processes  the  complex 
baseband  signal  Xn  =  xfnTg).  Like  x{f),  Xn  can  be  separated 
into  two  components: 


—  Xyy,l  T  Xxx^l 

Since  77„  and  are  real  processes,  their  autocorrelation 
functions  Vyy^i  and  rzz,i  are  both  real  and  even.  Thus,  rxx,i  is 
also  real  and  even.  ■ 

Lemma  2:  Let  a;„  be  a  lowpass  WSS  process  with  real  and 

even  autocorrelation  function  rxx,i  and  cutoff  frequency  Wc  G 
(0,7r/2).  Then,  there  exists  a  non-negative  integer  Nc  such 
that 

rxx,i  >  0  for  all|7|  <  (7) 

If  the  process  is  an  ideal  lowpass  process,  then  Nc  =  [tt/wj.]. 

Proof:  By  the  property  of  the  discrete-time  Fourier 
transform,  the  PSD  of  Xn  is  real  and  even.  Let  the  PSD  of 
Xn  be  defined  as 

{/(w),  0  <  w  <  Wc 

/(-w),  (8) 

0,  o.w. 

where  /(w)  is  a  positive  real  function.  Then, 

1 


^xx,l  — 


/(w)  COs{Lol)du! 


(9) 


where  s„  =  s{nTs),  Un  =  n{nTs),  and  wq  =  2'KfcTg.  The 
autocorrelation  function  of  s„  is  denoted  by  ;  =  rs{lTs). 
We  note  that  the  condition  {fb  <  fbw)  guarantees  that  s„  is 
non-white,  i.e.,  ^  where  Si  is  the  Kronecker  delta 

function. 

The  complex-valued  noise  component,  {un},  is  modeled 
as  a  circular  white  Gaussian  noise  process  with  mean  zero 
and  variance  Therefore,  the  autocorrelation  function  of 
the  noise  process  is  given  by  r„  ^  =  a'^Si.  Accordingly,  the 


We  observe  that  /(w)  cos(u;/)  >  0  for  all  lo  G  (0,7r/(27)). 
Thus,  Nc  =  [7r/(2wc)l  satisfies  (7).  Furthermore,  if  the 
process  is  an  ideal  “brick”  lowpass  process,  i.e.,  /(w)  =  1, 
then 


1 


^xxd  — 


cos{ujl)duj 


sin{u!cl) 

ttI 

Therefore,  Nc  =  [tt/wc]  satisfies  (7). 


(10) 


III.  Spectrum  Sensing  Algorithms 


Taking  N  samples  of  Xn,  i-e.,  x  =  ,xn-i),  a 

spectrum  sensing  algorithm  forms  a  decision  statistic  r(x) 
and  compares  it  to  threshold  A,  i.e., 


.  .  <  A  decide  Hq 
^  >  A  decide  Hi. 


(11) 


The  decision  statistics  in  [1]  as  well  as  our  proposed  detector 
are  based  on  the  estimates  of  autocorrelation  of  x. 

r  N-i-i 


A  )  7^  E  ^  0 

o  —  N  „=0 


(12) 


—r 


i ) 


I  <  0 


We  note  that  fi  is  an  unbiased  and  consistent  estimator  of 
Txx.i-  For  ease  of  notation,  in  the  following  we  will  drop  the 
dependence  of  the  decision  statistic  on  the  sample  data  x. 


A.  Zeng-Liang  signal  detector 

Zeng  and  Liang,  [1],  proposed  a  spectrum-sensing  technique 
base  on  the  following  decision  statistic. 


where  the  parameter  L  is  chosen  so  that  the  magnitude  of  the 
signal  autocorrelation  function  /|  is  significant  for  all  |^|  < 
L.  The  weighting  scheme  of  |f/p  on  the  left-hand  side  of  (13) 
is  formed  as  the  left-hand  side  is  derived  from  the  Frobenius 
norm  of  estimate  of  the  covariance  matrix.  Equivalently,  (13) 
can  be  reformulated  to  fit  (11)  as  follows. 


fzL  = 


(14) 


=  0  for  all  1.  On  the  other  hand,  under  Hq, 
=  0  for  all  I  0. 

In  general,  a  decision  statistic  is  designed  so  that  its  con¬ 
ditional  means  under  the  two  hypotheses  are  different.  The 
distinguishable  feature  in  t/  is  in  the  real  part  while  the 
imaginary  part  is  zero  under  both  hypotheses.  Hence,  we  form 
the  decision  statistic  for  our  proposed  detector  to  be 


Tin{uj)  =y^wi 
1=1 


To 


(17) 


Scaling  by  fg  results  in  a  constant  false-alarm  rate  (GEAR) 
detector.  The  limit  L  should  be  chosen  so  that  t;|^j  >0  for 
all  L  ^  Nc-  The  frequency  scanning  parameter  uj  enables  the 
algorithm  to  scan  across  the  frequency  band  in  order  to  locate 
the  center  frequency  of  the  captured  primary  signal.  While 
weighting  coefficients  wi  can  be  optimized,  in  this  paper  we 
have  opted  to  use  (15)  for  their  value. 


IV.  Performance  Analysis 

In  this  section,  the  performance  of  the  two  spectrum  sensing 
algorithms  presented  in  the  previous  section  is  evaluated 
analytically.  To  this  end,  we  first  determine  the  statistical 
distributions  of  the  autocorrelation  estimates  f/.  Then,  the 
performance  of  the  Zeng-Liang  detector  can  be  assessed 
following  the  procedures  used  by  Kay  [7]  to  analyze  the 
performance  of  an  autoregressive  detector.  The  performance 
of  our  proposed  detector  is  also  analyzed  using  the  same 
approach.  The  analysis  does  not  make  use  of  either  (Al)  or 
(A2)  and  only  assumes  that  {s„}  is  a  generic  complex-valued 
zero-mean  WSS  process  with  autocorrelation  function  rgj. 

A.  Asymptotic  conditional  distributions  of  fi 


with  weighting  function 


Wl 


A 


L  +  1-  |)| 
L  1 


(15) 


Because  the  autocorrelation  function  is  conjugate  symmetric 
and  the  (I  =  0)  term  is  always  1,  (without  affecting  the 
performance),  the  decision  statistic  in  (14)  can  be  simplified 
to 


^  A  \h\^ 

Tzl  = 


(16) 


B.  A  new  correlation-based  detector 


The  Zeng-Liang  detector  is  designed  to  account  for  the 
nonwhite  nature  of  the  (oversampled)  primary  communication 
signal.  By  further  incorporating  the  assumption  that  the  pri¬ 
mary  signal  Sn  is  (Al)  lowpass  and  (A2)  complex-valued  with 
independent  real  and  imaginary  components,  we  can  improve 
the  performance  of  the  autocorrelation-based  detector. 

By  (A2)  and  Lemma  1,  is  a  real- valued 

function.  Moreover,  introducing  (A2)  and  Lemma  2  reveal 
that  there  exists  an  integer  >  1  such  that  Xg  i  is  real 

and  strictly  positive  for  all  I  G  {—Nc,Nc).  Hence,  under 
Hi,  >  0  for  all  I  G  (-iV^,iVc)  while 


In  this  section,  we  present  the  probability  density  function 
(pdf)  of  the  autocorrelation  estimates  f/  in  (12).  Kay  [7]  has 
investigated  this  problem  for  the  case  of  an  autoregressive 
AR(1)  process  with  a  fixed  noise  power.  Here,  following  his 
approach,  we  present  a  generalization  to  his  results,  accounting 
for  an  arbitrary  SNR  value. 

To  compactly  formulate  the  statistics,  we  analyze  the  scaled 
version  of  f;,  namely 

yi  =  ^  (18) 

We  note  that  yi  can  replace  f;  in  (16)  and  (17)  without 
affecting  the  detection  statistics.  Let  yi  =  ai  jfdi.  Under 
each  hypothesis  {Hq  or  Hf)  and  for  large  N,  ai  and  /?/  can 
be  shown  to  be  jointly  Gaussian  distributed  by  the  central 
limit  theorem.  Consequently,  finding  the  conditional  mean  and 
conditional  (co)variances  of  oy  and  /3;  completely  determines 
their  statistics.  The  conditional  mean  under  each  hypothesis  is 
found  as  follows. 


E[ai\Hj^]  =  Si rjypr.i  (19) 

and 


E[Pi\Hr,]  =  r]yp^y 


(20) 


where  /rs,o  and  pi^i  =  9{rs,ie““'}/rs,o- 

Asymptotically  as  N  ^  oo  the  conditional  covariances  of  ai 
and  j3i  are  evaluated  as 

COv[Q;;,Q;m|-ffr,] 

1 

—  2^  H“  Pr,l-\-m)-)  (^l) 


cov[$i 

5  /5m 

1  ^ 

i^l  —  m  ^l^m)  “t”  {Pr,l—m  Pr,l+m)  5  (22) 


2Af 


'Af 


and 


COV  [A/ , I (.Pi^l  —  m  “t“  Pi^l+m)  •  (23) 


B.  Zeng-Liang  detector,  Tzl 

Zeng  and  Liang  [1]  only  provided  an  empirical  formulation 
for  the  probability  of  false  alarm  from  simulations  of  their 
detector  and  did  not  provide  the  probability  of  detection.  Using 
an  approach  similar  to  that  in  [7]  we  evaluate  the  cumulative 
distribution  function  (cdf)  of  the  decision  statistic  from  the 
characteristic  function  of  a  related  random  variable  [8]. 

First,  we  note  that  (16)  can  be  written  as 


T; 


ZL 


=  E 


Wl- 


1=1 


l«oP 


(24) 


Furthermore,  we  observe  that  the  cdf  of  can  be  written 
as 


Ftzl  (t) 


Let 

SzUt)  4  y^W(f)y  (26) 

where  ()^  denotes  matrix  transpose, 

y  =  [do  di  •  •  •  oiL  Pi  ■■■  Pl]'^  (27) 

and 

W(f)  =  diag([-f  wi  ...  WL  wi  ...  wl]).  (28) 

Then,  the  cdf  of  Tzl  can  be  expressed  in  terms  of  the  cdf  of 
SzL{t). 

FTzAt)=MSzL{t)<0} 

=  FszFtm 

where  FSzl  (i) ('^)  ts  the  cdf  of  Szl{^P  Hence,  the  cdf  of  Tzl 
can  be  evaluated  utilizing  the  characteristic  function  Pszl  A) 
of  SzL{t).  For  a  random  variable  X  with  characteristic 
function  (/)x(w),  its  cdf  can  be  computed  from  [8] 


Applying  this  characteristic  function  property  to  (29),  we  have 

FTzL{t)P)  =  FSzL{t)A) 

=  1  _  1  (31) 

2  TT  Jq  UJ 

Assuming  that  N  is  large,  op  and  Pi  are  jointly  Gaussian 
distributed  under  both  hypothesis.  Hence,  the  random  vector 
y  is  a  Gaussian  random  vector.  Denoting  y^  and  to  be 
the  conditional  mean  and  covariance  matrix,  respectively,  of 
y,  the  (conditional)  characteristic  function  of  a  quadratic  form, 
y^Wy,  of  Gaussian  vector  is  [9] 


Pz\Hr,A)  — 


|I-2ju;K^W|i/2 


X  exp 


jtuy^  W  (I-2jc.K,W)-'y, 


(32) 


We  note  that  there  is  no  known  closed-form  solution  to  the 
resulting  integral,  so  (31)  needs  to  be  numerically  evaluated. 

Finally,  we  can  formulate  the  probabilities  of  false  alarm 
and  detection  in  terms  of  Under  the  null  hypothesis 

Hq,  the  probability  of  false  alarm  given  the  threshold  Xzl  is 
computed  by 


Pfa.ZL  =  1  -  Ftzl\HoAzl) 

=  1  +  1  f°°  (33) 

2  TT  Jo 

Similarly,  under  alternate  hypothesis  Hi,  the  probability  of 
detection  given  the  threshold  Xzl  is  computed  by 


Pd,ZL  —  1  —  Ftzl\HiAzl) 

=  1  +  1  [°°  ^^i^SzLinPA} 

2  TT  Jq  UJ 


(34) 


Based  on  (19)-(23),  the  statistics  of  y  to  evaluate  (32)  can 
be  expressed  as  follows.  The  conditional  mean  of  y  is  found 
to  be 


y^=[l  0  •••  0]^ 

+  7[l  Pr-,1  •••  Pr,L  Pi,l  ■■■  (35) 

and  the  conditional  covariance  matrix  of  y  is  found  to  be 


K^=^diag([2  1 


Qi  Qt 
Q1  Qs 


•  (36) 


The  submatrices  Qi,  Q2,  and  Q3  in  (36)  are  composed  of  a 
Toeplitz  matrix  and  a  Hankel  matrix. 


Ql  —  2’(qr,0,L5  qr.0+)  + ‘^(qr,0+)  qr,L,2L),  (37) 


Q2  =  2'(qi^_i+_i,  qi^_i=L) 

+  P-Ai,!  ,L+l,^i,L+l,2L  ),  (38) 


Fx{x)  =  ^  - 

Z  TT 

^  /■“  lin[pxA)]  cos  LUX  -  Re[(()x(t^)]  sincua:^^^ 

Jo  ^ 


and 

Qs 


2’(qr,0,L-l5  qr,0,L-l) 

~  ^(qr,2,L-|-l)  qr+-|-l,2L),  (39) 


with 


Qr,a,6  —  Pr,a+1 

Pr,b\ 

(40) 

Pi,b\ 

(41) 

where  T(c,r)  represents  a  Toeplitz  matrix  with  the  first 
column  c  and  the  first  row  r,  7f(c,r)  represents  a  Hankel 
matrix  with  the  first  column  c  and  the  last  row  r.  Furthermore, 


and 


C.  The  proposed  detector,  Tjn 

To  follow  similar  steps  as  the  Zeng-Liang  algorithm,  we 
first  rewrite  the  proposed  decision  statistic  in  (17)  in  terms  of 
ai  and  jSi  by 

L 

Y,  wi[ai  cos(w()  +  /?/  sin(w()] 

Tjn{uj)  =  — - 


«0 


Let 


where 


SiN{t,uj)  =  w^(f)z(a;) 


(42) 

(43) 


z{uj)  =  [do  di  cos(w)  +  Pi  sin(u;)  •  •  • 

ctL  cos{coL)  +  Pl  sin(u;L)]^ 

and 

w(f)  =[-t  tui  •  •  •  wl\  ■ 

Then,  the  cdf  of  Tin{uj)  can  be  written  as 

=  Pr{S'/Ar(f,  w)  <  0} 

Similar  to  the  random  vector  y  in  (27),  the  random  vector 
z  is  Gaussian  under  both  hypothesis.  Hence,  S'/Ar(f)  in  (43) 
is  also  Gaussian  and  the  conditional  cdf  of  T/^v  is  given  by 

w'^(f)z^(w)  \ 


(44) 

(45) 

(46) 


^w'r(f)C^(w)w(f) 


(47) 


where  z^  is  the  conditional  mean  of  z,  C,,  is  the  conditional 
covariance  matrix  of  z,  and  Q{x)  is  the  Q-function.  Accord¬ 
ingly,  given  the  threshold  Xjn,  the  probability  of  false  alarm 
is  determined  to  be 


Pfa,IN{oj)  =  1  —  FTiN\Ho{^IN,i^) 

-w^(A7Ar)zoM 


=  Q 


(48) 


v/w^(A7Ar)Co(w)w(A/Ar)  J 

Similarly,  the  probability  of  detection  can  be  evaluated  by 

Pd,IN{i^)  =  1  — 

-w^(A7Ar)zi(w)  \  (49) 


=  Q 


■\/w^(A7Ar)Ci(w)w(A77v) 


The  conditional  statistics  of  z(w)  are  defined  as  follows. 
Under  H^j,  the  conditional  mean  z^(a;)  is  found  to  be 

z^(a;)=[l  0  •••  0]^ 

+  p[l  pi(w)  (50) 


P;(w)  =  Pr,i  cos(u;/)  +  p*,/  sin(a;0  (51) 

The  conditional  covariance  matrix  C^(w)  is  determined  by 

=  ^diag([2  1  •••  l]^)  +  77^ 

X  [T{po,l(w),Po.l(w)} +  H{po.l(w),Pl.2lH}]  (52) 
where 

Pa.bH  =  [Pa(w)  Pa+l(w)  •••  P&(w)]]^.  (53) 

The  expression  for  the  probability  of  false  alarm  in  (48)  can 
be  further  simplified  because  Cg  is  diagonal; 


The  false-alarm  rate  PfajN  is  independent  of  lu.  It  is  readily 
observed  from  (54)  that  the  probability  of  false  alarm  does  not 
depend  on  any  signal  parameters. 


V.  Numerical  Results 


In  this  section,  the  detector  performance  that  we  obtained  in 
the  previous  section  is  verified  against  Monte-Carlo  simulation 
results.  Following  common  configurations  are  used  throughout 
the  section.  All  detectors  use  the  same  number  of  samples, 
namely  N  =  1000.  The  detectors  are  configured  with  L  =  2. 
For  both  detectors  the  weighting  function  wi  is  as  defined  in 

(15). 

The  received  signal  contains  only  one  communication  sig¬ 
nal  under  Hi  case.  The  primary  signal  under  detection  is 
16QAM  signal  with  rectangular  pulse  shaping,  transmitted 
over  an  AWGN  channel.  The  transmitted  symbols  are  drawn 
randomly  with  equal  probabilities  among  all  possible  symbols. 
The  detectors  oversample  the  received  signal  at  Ns  =  3 
samples/symbol  (i.e.,  =  3/b).  Hence,  the  autocorrelation 

function  of  s„  is  given  by 


rs,i 


Ns-\l\ 

Ns  ’ 


0, 


\l\  <  Ns 

o.w. 


(55) 


In  the  simulations,  the  performances  are  evaluated  over  10000 
independent  trials  for  each  setup.  To  account  for  the  phase 
offset  6*0  in  (2),  6*o  is  randomly  drawn  from  [0,  27r)  for  each 
trial. 

Fig.  2  illustrates  the  receiver  operating  characteristic  (ROC) 
curve  of  the  two  detectors  under  fixed  7  =  — 12  dB,  and 
Fig.  3  shows  the  detection  probabilities  as  functions  of  the 
SNR  while  the  false-alarm  rate  is  fixed  to  a  constant  value 
of  0.01.  For  both  cases,  wg  =  0  and  Tjn  is  evaluated  with 
oj  =  ajQ  =  0.  The  figures  illustrate  the  improved  performance 
of  the  proposed  algorithm  over  that  in  [1].  Furthermore,  the 
theoretical  results  are  in  a  very  good  agreement  with  the  results 
from  simulation. 

Next,  we  observe  the  sensitivity  of  the  proposed  algorithm 
to  its  scanning  frequency  offset.  For  this  example,  we  fix 
Pfa  =  0.01  and  7  =  —8  dB  and  introduce  primary  carrier 


probability  of  false  alarm  Pf, 


Fig.  2.  ROC  curves  of  Tzl  and  Tjn(0}  (N  =  1000,  L  =  2,  7  =  -12 
dB,  Ns  =  3,  LOO  =  0). 


signal-to-noise  ratio  7  (dB) 


Fig.  3.  CFAR  Pd  vs.  7  (Pfa  =  0.01,  N  =  1000,  L  =  2,  Ns  =  3, 
loq  =  0). 

frequency  offset  ojq  =  0.1  tt.  Fig.  4  shows  the  probability 
of  detection  of  the  proposed  algorithm  as  a  function  of  the 
scanning  frequency  offset  uj.  The  Zeng-Liang  approach  which 
does  not  use  uj  is  also  shown  as  a  reference. 

As  designed,  the  detection  probability  Pd, in  peaks  at  uj  = 
ujQ  and  rolls  off  as  uj  moves  away  from  ujq.  The  detection 
rate  remains  high  over  the  signal  band  (roughly  2tt/Ns  rad 
around  wg),  indicating  that  the  proposed  method  can  identify 
the  spectral  location  of  s„.  This  is  an  advantage  over  the 
Zeng-Liang  algorithm  which,  by  construction,  cannot  detect 
the  location  of  s„  over  the  system  bandwidth. 


scan  frequency  w  ( x27r  rad./sample) 


Fig.  4.  CFAR  Pdjjv  vs.  u)  (Pfa  =  0.01,  N  =  1000,  L  =  2,  7  =  -8  dB, 
Ns  =  3,  OJQ  =  O.lV). 


VI.  Conclusion 

Autocorrelation-based  spectral  sensing  techniques  are  in¬ 
vestigated  and  analyzed.  A  new  technique  using  the  sample 
autocorrelation  function  of  the  received  signal  is  introduced 
and  its  performance  is  evaluated  theoretically.  Furthermore, 
the  analysis  is  extended  to  obtain  the  performance  of  the 
covariance-based  detector  in  [1].  The  performance  of  the  two 
detectors  are  also  obtained  through  Monte  Carlo  simulations. 
The  results  show  a  very  good  match  between  the  theoretical 
and  simulation  results.  This  verifies  the  accuracy  of  our 
asymptotic  analysis.  Furthermore  the  results  show  that  the 
proposed  method  outperforms  the  covariance-based  method  in 
[1]. 
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